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Abstract 



In recent work, it has been argued that multi-TeV masses for scalar su- 
perpartners are not unnatural. Indeed, they appear to have significant phe- 
nomenological virtues. Here we explore the implications of such 'focus point' 
supersymmetry for the dark matter problem. We find that constraints on 
relic densities do not place upper bounds on neutralino or scalar masses. We 
demonstrate that, in the specific context of minimal supergravity, a cosmo- 
logically stable mixed gaugino-Higgsino state emerges as an excellent, robust 
dark matter candidate. We estimate that, over a wide range of the unknown 
parameters, the spin- independent proton-neutralino cross sections fall in the 
range accessible to planned search experiments. 
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I. INTRODUCTION 



The idea that supersymmetry is an approximate symmetry of Nature, broken at a 'low' 
scale, is best and most concretely motivated by two general arguments. First, it addresses 
the question of why radiative corrections to the electroweak symmetry breaking scale, due 
to the exchange of virtual particles at some much larger unification or Planck scale, do 
not pull the electroweak scale to a much higher value (the gauge hierarchy problem Q). 
Second, it brings the unification of couplings calculation, which does not quite work if we 
use only the particles of the Standard Model, into adequate agreement with experiment [|| . 
However neither of these arguments, on the face of it, yields very specific constraints on the 
spectrum of masses of the supersymmetric (i?-parity odd) partners of ordinary (imparity 
even) particles: the first because it is inherently semi-quantitative, the second because it 
depends on these masses only logarithmically 

On the other hand, it is becoming notorious that the existence of superpartners with 
masses below 1 TeV is difficult to reconcile with established limits on various flavor violating 
effects, unless one assumes accurate degeneracy of squarks and of sleptons. Even if one 
assumes such degeneracy, supersymmetric contributions to CP violation and (in unified 
theories) proton decay become uncomfortably large, unless the squarks and sleptons are 
heavy. 

Therefore one is motivated to examine supersymmetric models with large scalar super- 
partner masses. The issue immediately arises: To what extent can such models, when 
critically scrutinized, embody the motivating features we mentioned in our opening para- 
graph? 

Recently it has been argued, using a seemingly reasonable objective definition of natu- 
ralness, that all squark and slepton masses can be taken well above 1 TeV with no loss of 
naturalness ||,f|] . The mathematical basis of this result is the existence of focus points in 
renormalization group trajectories, which render the weak scale (i.e., the Higgs potential) 
largely insensitive to variations in unknown supersymmetry parameters We should 

note that while in these models the squark and slepton masses are unusually large (com- 
pared to conventional wisdom), the electroweak gaugino and Higgsino particle masses are 
generically well below 1 TeV. 

In the focus point models, moreover, the quantitative aspect of the unification of coupling 
constants is actually slightly improved relative to more traditional supersymmetry models, 
reducing the need for significant high scale threshold corrections ||. 

Besides preserving and sharpening the original motivations for supersymmetry, the focus 
point models with heavy scalars ameliorate the problems noted above ||. Specifically, 
constraints on CP violation from electric dipole moments can be satisfied even with 0(1) 
phases in the fundamental parameters, and the predicted rate for nucleon decay — being 
proportional to the square of gaugino and inversely proportional to the fourth power of 
squark/slepton masses — is suppressed and less dangerous, though still perhaps accessible. 

In this note we will consider the implications of focus point supersymmetry for neutralino 
dark matter. We will concentrate on the focus point mechanism in its simplest 'default' 
incarnation, assuming a universal scalar mass within minimal supergravity. The focus point 
mechanism can also be realized in alternative models featuring gravity- 0,|J, gauge- [0], and 
anomaly-mediated || supersymmetry breaking. In addition, the requirement of a universal 
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scalar mass, though sufficient, is not necessary. For the precise set of requirements, see 
Ref. 1. 

Of course, the cosmology of minimal supergravity (with universal scalar mass m < 
1 TeV) has been studied extensively. Prominent among conclusions one finds in the literature 
are the following: 

(i) In essentially all of parameter space, the lightest neutralino \ is Bino-like, with gaugino 
fraction R x > 0.9 |5|. (See below for the definition of R x .) 

(ii) The neutralino relic abundance Q x h 2 increases as mo increases. The requirement that 
neutralinos not overdose the universe, along with the requirement that the lightest super- 
symmetric particle (LSP) be neutral, then typically leads to stringent upper bounds on 
scalar and/or gaugino masses, independent of naturalness considerations []T0|— p79[] . 

(iii) Dark matter detection rates decrease as mo increases. The rates predicted range widely 
over parameter space, but the prospects for discovery at ongoing experiments grow increas- 
ingly dim as m increases 



In our analysis, we confirm these conclusions — for m < 1 TeV. However, we find that 
they fail for focus point models with m > 1 TeV. Instead, we find that in such models: 

(i) The lightest neutralino is a gaugino-Higgsino mixture over much of parameter space. 

(ii) Neutralino relic abundances do not overdose the universe. Rather, they tend to lie in 
the cosmologically interesting range. Therefore, in particular, cosmological considerations 
alone do not place stringent upper limits on superpartner masses, even apart from possible 
conspiracies involving poles |21|J2^] or co-annihilation ||l6p3" |. 



(iii) Direct detection rates may be fairly large, with proton-neutralino cross sections plausibly 
falling within the range 10 -6 -10 -8 pb. (See the more detailed discussion and plots below.) 

These results are especially interesting since ongoing and planned dark matter detection 
experiments promise to cover the indicated range. 



II. NEUTRALINOS IN FOCUS POINT SUPERSYMMETRY 

The models we consider, like many supersymmetric models, contain a natural cold dark 
matter candidate P^j . Assuming that R = (— S ) B + L + 2S is accurately conserved, the lightest 
supersymmetric (.R-odd) particle is stable. We will be considering models where, with no 
special adjustment, this is a neutral weakly interacting massive particle (WIMP), which is 
commonly called the neutralino. 

The relic densities and detection rates of neutralino dark matter can be calculated, assum- 



ing straightforward extrapolation of Big Bang cosmology to T ~ 10 GeV ||25|| . In principle 
they depend on details of the entire supersymmetric spectrum, but in practice, especially 
in focus point models, they are mainly determined by the properties of the neutralino LSP 
itself. 

Assuming i?-parity conservation and minimal field content, at tree level the neutralino 
mass matrix is 
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/ Mi — mz cos (3 sin 9w mz sin (3 sin 9w \ 

M 2 m z cos j3 cos 9w — mz sin j3 cos 9w 

—mz cos (3 sin mz cos (3 cos 9w — [A 

V mz sin P sm9w — mz sin /3 cos 9w / 



(1) 



in the basis (—iB,—iW 3 ,H®,H®). Here Mi and M2 are the soft Bino and Wino masses, 
/i is the Higgsino mass parameter, and tan (3 = (H®) / (H%) is the ratio of Higgs vacuum 
expectation values. We parametrize the gaugino/ Higgsino content of the lightest neutralino 
according to 

X = ai (-iB) + a 2 (-iW 3 ) + a 3 H° + a 4 H° d , (2) 

and define the gaugino fraction of x to be 

R x = | a,! | 2 + | a,2 1 2 • (3) 

In the framework of minimal supergravity, there are 4 continuous parameters and 1 
binary choice: 

m , Mi/2, A , tan/3, sgnQi) . (4) 

Here, m , Mi/ 2 , and A Q are, respectively, the universal scalar mass, gaugino mass, and 
trilinear scalar coupling at the grand unified theory (GUT) scale Mgut — 2 x 10 16 GeV. 
Given values for these input parameters, all the couplings and masses of the weak scale 
Lagrangian are determined through renormalization group (RG) evolution. In our work 



we use two-loop RG equations |26[] and include one-loop threshold corrections from super- 
symmetric particles to the gauge and Yukawa coupling constants [^,^]. We minimize the 
Higgs potential after including all one-loop effects, and include one-loop corrections in all 



superpartner masses ||28|| . This framework imposes several specific relations among the weak 
scale supersymmetry parameters. In particular, Mi ~ M 2 /2 ~ 0.4Mi/ 2 , and is fixed by 
the condition of electroweak symmetry breaking, which, at tree-level, is 

1 2 m\ -m\ u tan 2 /? 222 (< s 
-m 7 = 7T- 7^ u ~ — m H — a , (5) 

2 z tan 2 (3-1 p Hu p ' v ' 

where the last relation is valid for moderate and large tan (3. 

The essence of the focus point mechanism is the observation that the weak scale value of 
m 2 H is remarkably insensitive to variations in the fundamental GUT scale supersymmetry 
parameters, even for multi-TeV m . By Eq. (||), therefore, for tan/3 > 5 [|J, the electroweak 
scale is insensitive to variations in these parameters, and in this sense multi-TeV values of 
mo are natural. 

Contours of gaugino fraction R x are presented in Fig. [I]. (For reference, contours for 
the fine-tuning parameter c, as defined in Ref. 0], are also given in Fig. |l|. While large 
M1/2 leads to large fine-tuning, large mo does not, as a result of the focus point discussed 
above.) For m < 1 TeV, the lightest neutralino is nearly pure Bino, with R x > 0.9. This 
well-known result arises from the circumstance that RG evolution typically drives m 2 H large 
(relative to the gaugino masses) and negative, which by Eq. (|) implies much larger than 
Mi and M 2 . 
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FIG. 1. Contours of constant LSP gaugino fraction R x (in percent) in the (mo, M]y 2 ) plane for 
A = 0, fi > 0, and two representative values of tan/3. The shaded regions are excluded by the 
requirement that the LSP be neutral (top left) and by the chargino mass limit of 95 GeV (bottom 
and right). Dashed contours are for the fine tuning parameter c = 20, 30, 50, 75 and 100, from 
below (see text). 



For m > 1 TeV, however, \ contains significant amounts of Higgsino over much of the 
parameter space. For large and moderate values of tan /3, large values of m generate positive 
corrections to m 2 Hu , rendering it less negative. Thus if all other parameters are fixed, we find 
that as m increases, decreases, and eventually \p\ ~ M 1; M 2 , which leads to significant 
mixing between Higgsino and gaugino states. As mo increases further, the Higgsino content 
of x increases until ultimately one enters the domain < 95 GeV, which is excluded by 



limits from chargino searches at LEP II [29]. At that point, we have entered the shaded 
region of Fig. [I]. 

Note that the above discussion holds for tan (3 > 5 where the focus point mechanism 
is operative. For small tan/3, p becomes sensitive to m? H , and large mo gives negative 
contributions to m% . Both of these effects imply that, as mo increases, \p\ also increases, 
and so there is no mixed gaugino-Higgsino region for small tan (3. 



III. RELIC ABUNDANCE 

There is ample and increasingly precise evidence that the energy density of the universe is 



not dominated by the observed luminous matter ||30|| . The evidence is conveniently expressed 
in terms of ratios fli = Pi/p c , where the pi are energy densities and p c is the critical density, 
and h « 0.65±0.1, the Hubble parameter in units of 100 km s _1 Mpc -1 . The luminous matter 
density is roughly fl\ um ~ 0.005, and the successful predictions of Big Bang nucleosynthesis 
require baryon density Q^h 2 fa 0.02. At the same time, rotation curves of spiral galaxies 
require matter density fl m ^ 0.1, and a variety of observations, ranging from the velocities 
of galaxies within galactic clusters to the luminosities of Type la supernovae, favor values 
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FIG. 2. Contours of constant relic density £l x h? in the (mo,M^ 2 ) plane for Aq = 0, \i > 0, 
and two representative values of tan/?. In the black shaded region, \2m x — < 5 GeV, and the 
h pole becomes important. The light shaded regions are as in Fig. |l[ 



in the range 0.2 < Q m < 0.6. There is therefore a need for both baryonic and non-baryonic 
dark matter. For non-baryonic dark matter, which is of interest here, galactic rotation curves 
require Odm/i 2 > 0.025, and large scale measurements suggest 0.1 < ^dm^ 2 ~ 0.3. 

Neutralinos remain in thermal equilibrium until their pair annihilation rate drops below 
the Hubble expansion rate. The relic density of neutralinos is therefore determined primarily 
by their cross section for pair annihilation in the early universe. Neutralinos annihilate to 
fermion pairs through the t-channel exchange of sfermions /, and via XX —> Z, A, h, H —>■ ff, 
where A is the CP-odd Higgs boson, and h and H are the lighter and heavier CP-even Higgs 
bosons. When kinematically accessible, annihilation to two body states of gauge and Higgs 
bosons are also relevant. 

The calculation of neutralino relic density has been discussed at length in the literature. 
We use the results of Drees and Nojiri PT| , |3~T| , as implemented by Jungman, Kamionkowski, 
and Griest 
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S- and P-wave contributions to all tree-level processes with two-body final 
states are included. 

In Fig. we display the neutralino relic density fi x h 2 in the (m , Mi/ 2 ) plane for two 
representative values of tan/?, A = 0, and fi > 0. These plots exhibit the most important 
parameter dependencies. Variations in Aq within its natural range have little effect, as the 
A parameter flows to a fixed point at the weak scale. Similarly, changing the sign of /i does 
not produce qualitatively new results. We have chosen the sign less constrained by b — > sj 
(see below). 

In most of parameter space, the expected accuracy in Fig. |2| is O(10%) |32|]. More 
significant errors are possible in special regions. First, a more refined treatment is necessary 
just below annihilation thresholds P3|J3^ ]. Such a treatment will smooth out the kinks in 
the contours of Fig. ^, which are caused by the opening of annihilation channels such as 
WW and ZZ. Second, co-annihilation is important very near the left and right borders of 



the allowed region, where the LSP is nearly degenerate with staus |T6|,^5j and Higgsinos [|36 
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respectively. Finally, s-channel poles also require a more careful analysis [p3|,p4 | p7 , p^ . The 
Z pole possibility, 2m x = mz, is now essentially eliminated by the chargino mass bound of 
95 GeV p9| , but a h pole is possible near the bottom of the allowed region, as indicated 
in Fig. H and A poles, though possibly significant for high tan f3 in the m < 1 TeV 
region [p2"p8f , are absent for mo > 1 TeV for all tan/? < 50 and fi > 0. (The sign of 
/i enters through the finite supersymmetric threshold corrections to the bottom Yukawa 
coupling PjJ.) Thus, although special effects warrant more sophisticated analysis in limited 
regions of parameter space, they are absent in the bulk of the focus point region. They do 
not affect our main conclusions. 

For sub-TeV m , we see that Q x h 2 is a monotonically increasing function of m . In 
this region, the LSP is nearly Bino-like, and so the dominant process is XX ~^ ff through 
i-channel / exchange. As mo and mj increase, this annihilation process is suppressed, 
and the relic density grows. This behavior, along with the requirement of a neutral LSP, 
has led many authors to conclude that upper bounds on the relic density typically impose 
stringent upper limits on superpartner masses pOHTTj. As these results are independent of, 
and less subjective than, considerations based on naturalness or the fine structure of gauge 
coupling unification, they have often been taken, by the authors themselves or others, as 
robust phenomenological upper bounds. If correct, this conclusion would have important 
implications for supersymmetry searches at the Tevatron and the LHC, and in planning for 
future linear colliders JI^,|18|,|19| . 

But when we consider mo > 1 TeV, motivated by the focus point or otherwise, we find 
that the behavior of this function reverses. The reason for this behavior is not difficult 
to locate. For although as mo increases the t-channel sfermion exchange process is more 
and more suppressed, as noted in Sec. [H|, at the same time the LSP gradually acquires a 
significant Higgsino component. Because of this, other diagrams become unsuppressed, and 
pair annihilation again becomes efficient. The cosmologically interesting region, with a mixed 
gaugino- Higgsino LSP and 0.1 ^5 ^\^ 2 ~ 0-3; extends at least to mo ~ 10 TeV, Mi/ 2 ~ 
6 TeV and m x ~ 2.5 TeV. Thus, for all collider applications, cosmological considerations 
do not provide upper bounds on superpartner masses, even in the constrained framework 
of minimal supergravity. It has been noted previously that cosmological bounds might be 
weakened or eliminated by the effects of co-annihilation HT6||23ll or poles pl|,22|. Here, we see 
that even without appealing to such conspiratorial effects, sweeping claims that cosmology 
provides stringent upper bounds on superpartner masses are unfounded. Note that the 
cosmologically interesting region of parameter space with m > 1 TeV is in no sense small. 
For fixed M 1 / 2 , this region extends over one to several hundred GeV in m , depending on 
tan /3. It is comparable in area to the allowed region with m < 1 TeV. 

For all tan/5 > 5, we find relic densities in the cosmologically interesting range 0.025 < 
Q x h 2 < 0.3 over most of the focus point regime with mo > 1.5 TeV. In fact, as can be seen 
from Fig. |2|, for tan/5 = 50, Q x h 2 never saturates the upper bound, even for m ~ 1 TeV. 
This additional suppression of Q x h 2 is not due to any enhancement of the LSP's Higgsino 
component from large tan/3. As can be seen in Fig. [I], the gaugino purity of the LSP 
is largely independent of tan/3. Rather it results from two enhancements to the process 
XX ~~ * A — > //. First, the Aff coupling is proportional to tan/3 for isospin — | fermions. 
Second, for large tan/3, the bottom Yukawa coupling hf, becomes significant and lowers m^. 
This can be seen by noting that, for h t = hf,, symmetry under interchange t <-> b implies 
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that m 2 Hd also has a weak scale focus point, and so htla ~ 0(100 GeV) <C m . 

For small tan j3, as noted in Sec. 0, there is no mixed gaugino-Higgsino region for large 
m . Thus, fl x h 2 grows monotonically as m increases. Low values of tan (3, such as tan (3 = 2, 
have been widely considered in the literature. They are, however, increasingly disfavored by 
current null results in Higgs searches, and, in any case, are less natural than moderate and 
large values, as low tan/3 typically leads to large values of \n\ f|. 



IV. DIRECT DETECTION 

Having shown that focus point models predict cosmologically interesting densities of dark 
matter, we now consider their predictions for dark matter detection rates. The interactions 



of neutralinos with matter are dominated by scalar (spin-independent) couplings pl| , |32 
These interactions are mediated either by t-channel h and H exchange, or via %q — > q — > XQ- 
The former diagrams contain h\X an d H\X vertices, which are suppressed for Bino-like x- 
They also contain h, if -nucleon vertices. These are dominated by contributions arising from 
fundamental couplings of h, H to the strange quark and indirectly to gluons through heavy 
quark loops. 

In the calculation of dark matter detection rates, we assume a Maxwellian x velocity 
distribution with velocity dispersion v = 270 km/s, and a local dark matter density of 
po = 0.3 GeV/cm 3 . 

A major uncertainty in the estimate of dark matter detection rates arises from the poorly 
determined value of the nucleon matrix element 

f T . = (N\m s ss\N)/m N . (6) 

If it is correct to treat the effect of the strange quark mass term m s ss in first order pertur- 
bation theory, then we can identify Jt s as the fractional change in the nucleon mass from 
what it would be in a world with massless s quarks. Values as high as Jt s = 0.62 and as low 
as fx 3 = 0.08 have been considered in the literature [ 40fl . 



The traditional strategy for evaluating fx a is to combine an 'experimental' determina- 
tion of the cr-term (N\ m "+ md (uu + dd)\N) with a fit of the flavor S77(3)-breaking term 
(N\ |m s (2ss — uu — dd) | N) to the baryon octet splitting, and finally the ratios of light quark 
masses from the n/K mass ratio. This procedure is problematic, however, in several ways. 
First, theory relates the value of the a term to the value of the irN scattering amplitude 
at an unphysical point, and a rather elaborate and numerically unstable extrapolation is 
necessary (hence the quotation marks). Second, in getting to the quantity of interest, the 
uncertainty in the a term is amplified by the large factor 2m s / (m„ + m^). Third, the nucleon 
matrix element of the 577(3) breaking term is presumably dominated by the valence quark 
terms, so that the term of interest is subdominant. Fourth, and exacerbating all these, the 
use of first-order perturbation theory for the baryon splittings is questionable to begin with. 

Clearly it would be very desirable to have a first-principles estimate of Jt s - Straightfor- 
ward evaluation in lattice gauge theory appears very difficult, because three-point amplitudes 
with disconnected components are noisy and computationally demanding. Here we would 
like to suggest an alternative strategy that avoids any measurement of three-point ampli- 
tudes. The basic idea is to vary the strange quark mass, starting with a very large value and 
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bringing it down by small increments. In the heavy quark limit we have a universal answer 
by integrating out the heavy quark to obtain a gluon operator, which is (essentially) the 
trace of the energy- momentum tensor with a definite coefficient |41]|. And if we change the 



strange quark mass by a small amount Am s , we have for the induced change in the nucleon 
mass 

Am N = (N\Am s ss\N) , (7) 

by first-order perturbation theory. Of course, it is implicit in this procedure that one is doing 
a fully unquenched calculation, with dynamical quarks (including the strange quark). Since 
the error is at least quadratic in the step size, one should achieve good accuracy by taking 
the steps sufficiently small. Note that, since the nucleon state is changing at each step of 
the calculation, it is not correct to do the whole thing in one jump, despite the apparent 
linearity of Eq. (|7|). In any case, Eq. (|7|) allows one to extract the matrix element of interest 
by measuring masses only, without ever involving three-point functions. 

No result of this type is currently available, and so at present one must rely on informed 
guesswork for a numerical estimate of fx s ■ There is a very extensive literature on the subject, 
which we will not review here. We will confine ourselves to two brief observations. 

If one treats the strange quark as heavy, rather than light, then one obtains the universal 
value /r Hcav = 0.074. This is presumably a lower limit. But of course if one takes the strange 
quark mass to zero, then Jt s will likewise vanish. Thus, in considering values of fx s that are 
many times the heavy quark value, one is implicitly postulating that this quantity varies 
very rapidly as a function of m s , which may be implausible. 

Quenched calculations of the nucleon to rho mass ratio, and even more so calculations 
with two dynamical quarks, agree remarkably well with experiment, down to the few percent 
level. This fact can only be consistent with the idea that the strange quark is significantly 
perturbing the nucleon mass if the strange quark contributes an accurately equal proportion 
of the rho mass. In general, that would seem to require a conspiracy. Just such a "conspir- 
acy" does take place for heavy quarks, since they can be integrated out in favor of a charge 
renormalization, whose primary effect is just to change the overall mass scale (dimensional 
transmutation). However, as we have seen, if it is valid to treat the strange quark as heavy, 
then its contribution is small. 

While these considerations are far from definitive, they do seem to us to make the larger 
end of the range of values discussed for f Ts appear dubious. In preparing our plots, we 
have adopted the conventional value /y a = 0.14. For extreme values in the range quoted 
above, the proton-neutralino cross sections plotted may decrease by roughly a factor of 2, 
or increase by a factor of 5. 

In Fig. ^, we present proton-neutralino cross sections <jp. For m ^ 1 TeV, the dominant 
contribution is through squarks, and so the interaction rate decreases monotonically for 



increasing m [20]. A naive extrapolation would then suggest very small detection rates in 
the focus point region. However, for mo ^ 1 TeV, this behavior is reversed — in this region, 
as mo increases, the LSP's Higgsino content increases, and so the Higgs boson diagrams 
become less and less suppressed. (This feature was also noticed in Ref. W^i where the 
region of large mo was singled out on the basis of Yukawa coupling unification in SO(10) GUT 
models.) Detection rates may therefore be large in the focus point region. This is especially 



true for large tan/3, where the effect noted in Sec. Ill, a relatively light H boson with tan/3 
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FIG. 3. Contours of proton- x cross section ap in pb in the (mo, Mi/ 2 ) plane for Aq 
and two representative values of tan (3. The shaded regions are as in Fig. [|. 



0, \i > 0, 



enhanced Hff couplings, leads to large and dominant heavy Higgs amplitudes ||3 . 

In Fig. |], we show the correlation between relic density and proton-neutralino cross 
section for a representative sample of points in parameter space. As usual, there is a strong 
anti-correlation, with large o> corresponding to low Q x h 2 . Nevertheless, we find that, for a 
given tan (3 and Q x h 2 , points with large mo have proton-neutralino cross sections comparable 
to those with conventional sub-TeV mo- For the optimal case of large tan/?, points with 
tuq > 1 TeV and Q x h 2 > 0.1 (Q x h 2 > 0.025) have cross sections as large as ap ~ 2 x 10~ 7 
pb (a P ~ 10~ 6 pb). 



V. DISCUSSION 

Several ongoing and planned experiments, using a variety of techniques, are devoted to 
the search for WIMP dark matter |P|]. Assuming that the matter- WIMP interactions are 
dominated by spin-independent couplings, these experiments may be compared by scaling 
the matter- WIMP cross sections for each detector material to the proton- WIMP cross section 
ap and displaying the sensitivity in the (m x , ap) plane. 

In Fig. |5| we plot predictions of focus point supersymmetry in the (m x ,ap) plane. We 
find that, for a given m x , cross sections are typically maximized for large m , and large 
cross sections near current sensitivities are possible. At present, the DAMA Collaboration 
has reported data favoring the existence of a WIMP signal in their search for annual mod- 



ulation When WIMP velocity uncertainties are included [I6[], the preferred range of 



parameters is 10" 6 pb < a P < 10~ 5 pb and 30 GeV < m x < 200 GeV (3a CL). The DAMA 



result has been criticized ||47[| , and a recent analysis of the CDMS Collaboration excludes 



the 3a DAMA region at > 84% CL [jiS ]. While the situation is at present unclear, in the 
near future both experiments will improve their sensitivities significantly. It is noteworthy, 
in any case, that current experiments are on the verge of probing the region of parameter 
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FIG. 4. Regions of the (Q, x h 2 , ap) plane populated by minimal supergravity models with 
rriQ < 1 TeV (yellow, light), 1 TeV < ttiq < 1.5 TeV (green, medium), and 1.5 TeV < too (blue, 
dark). The parameters scanned are those given in Fig. |l], subject to the additional naturalness 
constraint Mi/ 2 < 400 GeV. We assume neutralino velocity dispersion v = 270 km/s and local 
density po = 0.3 GeV/cm 3 , and fx s = 0.14 (see text). 



space arising in focus point models. 

It is also interesting to compare the discovery potentials of dark matter searches and 
future colliders in focus point models. In the focus point scenario, all squarks and sleptons 
are beyond the reach of the Tevatron and may present significant challenges even for the 
LHC. However, the lightest Higgs boson has mass m h < 120 GeV, and so may be discovered 
in Run II of the Tevatron. The correlation between and op is presented in Fig. [| We 
note that focus point models provide a natural setting for Higgs masses above the current 
bound, as the large squark masses raise the Higgs mass through their radiative corrections. 
At the same time, both Higgs and dark matter searches are promising, and will confront the 
focus point scenario over the next few years. 

It has been argued that the regions of parameter space which give large dark matter 



detection rates may also predict large rates for other processes, such as proton decay |49 
and b — > 57 |50[. Conversely, present bounds from such processes may exclude the largest 
possible rates for dark matter detection. The proton decay rate, as noted in the introduction, 
scales as m^ 4 , so it is relatively suppressed in the focus point models. The process b — ► 57 
provides a stringent constraint on large tan/? only for one sign of \i. For the positive sign 
(the one adopted in our plots), destructive interference between the chargino and charged 
Higgs diagrams typically leads to predictions within the experimental bounds. 

Focus point scenarios can also imply interesting rates for anti-matter detectors and ex- 
periments sensitive to dark matter annihilation in the cores of the earth or sun E3| . 



Finally, of course, neutralinos and charginos themselves are targets for discovery at collid- 
ers. If they are found, detailed measurements of their masses and compositions could then be 
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made, and might even be used to exclude neutralinos as dark matter candidates ||51|| . How- 
ever, colliders will never be able to distinguish unstable particles with long decay lengths 
from those stable on cosmological time scales. Since some well-motivated models with low 
energy supersymmetry breaking predict LSPs which are stable on accelerator detector, but 
unstable on cosmological, scales, in principle identification of neutralinos as the dark mat- 
ter of our universe cannot be made at accelerators. Ultimately this identification requires 
detection, either directly or indirectly, of neutralinos permeating space. 
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